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We present a method to obtain higher order integrals and polynomial al- 
gebras for two-dimensional superintegrable systems from creation and anni- 
hilation operators. All potentials with a second and a third order integrals 
of motion separable in Cartesian coordinates were studied. The integrals of 
motion of two of them do not generate a cubic algebra. We construct for 
these Hamiltonians a higher order polynomial algebra from the creation and 
annihilation operators. We obtain quintic and seventh order polynomial alge- 
bras. We give also for polynomial algebras of order 7 realizations in terms of 
deformed oscillator algebras. These realizations and finite dimensional uni- 
tary representations allow us to obtain the energy spectrum. 



1 Introduction 

Over the years many articles were devoted to superintegrability [1-12]. In 
classical mechanics a Hamiltonian system with Hamiltonian H and integrals 
of motion Xa 

H = ]^gikPiPk + V{x,p), Xa = fa{x,p), a=l,...,n-l , (1.1) 

is called completely integrable (or Liouville integrable) if it allows n inte- 
grals of motion (including the Hamiltonian) that are well defined functions 
on phase space, are in involution {H,Xa}p = 0, {Xa,Xb}p = 0, a,b=l,...,n-l 
and are functionally independent ({, }p is a Poisson bracket). A system is su- 
perintegrable if it is integrable and allows further integrals of motion Yh{x,p), 
{H,Yb}p = 0, b=l,...,k that are also well defined functions on phase space 
and the integrals{i/, Xi, Xn-i, Yi, Yk} are functionally independent. A 
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system is maximally superintegrable if the set contains 2n-l functions such 
integrals. The integrals Yf, are not required to be in evolution with Xi,...Xn-i, 
nor with each other. The same definitions apply in quantum mechanics but 
{H, Xa, Yf,} are well defined quantum mechanical operators, assumed to form 
an algebraically independent set. Superintegrable systems in classical and 
quantum mechanics possess many properties. These systems appear to be 
important from the point of view of mathematics and physics. One of their 
interesting property are their nonabelian algebraic structure generated by 
their integrals of motion. This algebraic structure can be a finite dimensio- 
nal Lie algebra [1,2,3], a Kac-Moody algebra [13] or a polynomial algebra 
[14,15,16]. These polynomial algebras were related to deformed oscillator al- 
gebras and parafermionic algebras. Superintegrable systems are also related 
to systems studied in supersymmetric quantum mechanics [17-25]. Super- 
symmetry in quantum mechanics allows us to obtain the wave functions and 
the energy spectrum. 

This article follows a series of articles [26,27,28,29,30,31] devoted to super- 
integrable systems in classical and quantum mechanics with third order in- 
tegrals. In two-dimensional Euclidean space E2 there are fourteen quantum 
systems with a second and a third order integrals [27]. These systems were 
investigated from the point of view of cubic algebras and supersymmetric 
quantum mechanics [30,31]. The Hamiltonians and the integrals of motion 
were of the following form 

H=^ + !f + g,{x)+g,{y) , (1.2) 

A=^-!f + 9i{x)-g2{y) , (1.3) 
B= J2 Aijk{Llp{p'^} + {h{x,y),pi} + {l2{x,y),p2} , (1.4) 

where {, } is an anticommutator, L3 = XP2 — yPi is the angular momentum 
and [H,A] = [H,B]=0. The constants Aij^ and functions li and I2 are known 
[27]. 

We considered the most general cubic algebra generated by the integrals 

[A,B] = C (1.5a) 
[A,C] = aA'^ + /3{A,B} + -fA + 6B + € (1.5b) 
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[B,C]=ijA^ + uA'^-f3B^-a{A,B} + ^A--fB + C ■ (1.5c) 

where a, /?, 7, S, e, ^, v and C are polynomials in H. We constructed realiza- 
tions in terms of deformed oscillator algebras and found Fock type represen- 
tations [29,30]. For twelve cases, the cubic algebra belongs in the following 
particular case 

\A,B\=C, \A,C\ = ^B, \B,C\= iiA^ ^vA^ ^iA^C, . (1.6) 

There are 6 irreducible (i.e. the third order integral is not a consequence of 
lower order one) quantum superintegrable Hamiltonians with a second and 
a third order integrals separable in Cartesian coordinates written with a ra- 



tional function : 


Potential 1. 


V = 


Potential 2. 


V = 


Potential 3. 


V = 


Potential 4. 


V = 


Potential 5. 


V = 


Potential 6. 


V = 



f^2\9x^+y^ I 1 I 1 1 
\-~8^ ^ (s/-a)2 (j/+a)2J 

^ (8^[(^ + ^ ) F ~^ (a;+a)2 + (x-a)2] 



(j/+a)2 (y—a)^ (x+a)^ {x—a)^J 

The integrals A,B,C respectively of order 2,3 and 4 of the Potential 5 and 
6 do not generate a cubic algebra [30]. We studied these six potentials from 
the point of view of supersymmetric quantum mechanics (SUSYQM) and 
obtained the ladder operators, wave functions and energy spectrum. 
Let us present the organization of this paper. In Section 2, we will show 
how we can generate higher order integrals and a polynomial algebras for 
two-dimensional Hamiltonians constructed from two one-dimensional Hamil- 
tonians and their creation and annihilation operators. We present the poly- 
nomial algebra for two cases. In Section 3, we apply the results of Section 
2 to the Smorodinsky-Winternitz potential and the Potential 5 and 6. The 
method presented in Section 2 can also be used to generate new superinte- 
grable systems from known one-dimensional Hamiltonians for which creation 
and annihilation operators polynomial in momenta exist. We generalize the 
potential with a fourth Painleve transcendent given by Eq.(l.lO) and also 
construct its polynomial algebra with the results of Section 2. These results 
extend the number of known superintegrable systems involving the Painleve 
transcendents. In Section 4, we give for a class of polynomial algebras of or- 
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der 7 the realizations in terms of deformed oscillator algebras. In Section 5, 
we use the results of Section 3 and 4 to obtain the Fock type unitary repre- 
sentations and the corresponding energy spectrum of the Potential 5 and 6. 



2 Polynomial algebras 

Let us consider a two-dimensional Hamiltonian separable in Cartesian 
coordinates 

H{x, y, P^, Py) = Hi{x, P.,) + H2{y, Py), (2.1) 

for which creation and annihilation operators A^, Al, Ay and A^y (polyno- 
mials in momenta) exist. These operators satisfy the relations 

[H„Al] = KAl [H,,Al] = XyAl . (2.2) 
The following operators 

h = Al^A-y, f2 = A^Al- , (2.3) 
commute with the Hamiltonian H given by Eq.(2.1) 

[H,h] = [H,f2]=0, (2.4) 

if 

mX^ — nXy = 0, m, n G Z"*" . (2.5) 

Creation and annihilation operators allow to construct polynomial integrals 
of motion. The following sums are also polynomial integrals that commute 
with the Hamiltonian H 

h = Al^A;-A^Al\ h = Al-A; + A:^Al- . (2.6) 

The order of these integrals of motion depends of the order of the creation 
and annihilation operators. The separation of variable in Cartesian coordi- 
nates implies the existence of a second order integral K = — Hy. The 
creation and annihilation operators can also provide a method to determine 
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the polynomial algebra of the two-dimensional superintegrable systems gi- 
ven by Eq.(2.1). We require that operators and Al satisfy these further 
relations 

[A,,Al] = P{H,) = Q{H, + A,) - Q{H,), (2.7) 

[Ay,Al] = R{Hy) = S{Hy + A,) - S{Hy), (2.8) 

where P{Hx) is of order Ux and Q{Hx) is of order rix + 1 and R{Hy) is of 
order Uy and S{Hy) is of order Uy + 1. The algebras given by Eq.(2.7) and 
(2.8) are deformed oscillator algebras. They can be interpreted polynomial 
superalgebras i.e. {A^, A^} = Q{Hx + Xx) + Q{Hx). 

2.1 Case A^; = Ay = A 

We consider the case 

\x = \y = A, (2.9) 
and we take the following linear combination 

A = 2{Hx - Hy), h = {AlAy - AxAl), h = 4A(4.A, + A.Aj). (2.10) 

We use the integrals given by Eq.(2.10) to construct polynomial algebras of 
the Hamiltonian given by Eq(2.1). 

[A/i]=/2, [A/2] = 16A2/i, (2.11) 
[/i, h] = 8A(g(i(iJ + ^Am^-iH - ^A) + A) 

-Qi^iH +^-A) + X)Si^iH -^A)). 

The order of this polynomial algebra depends of the order of the polynomial 
Q and S. An example of such construction was used write the angular momen- 
tum algebra as two independent harmonic oscillators [32]. The construction 
of integrals of motion from creation and annihilation operators was discus- 
sed in earlier articles for specific examples : The two-dimensional harmonic 
oscillator [33], the three dimensional harmonic oscillator and the hydrogen 
atom [34], the Smorodinsky-Winternitz potentials [6,35] and the anisitropic 
oscillator in Ref.36. 
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A, ^ 1 

Xy 2 

We consider the case 



2.2 CclSG — r) 

Ay I 



2X., = Xy = X , (2.12) 
and take the following integrals 

A = 2{H^ - Hy), h = {Al'Ay - AlAl),h = ^XiAl'Ay + A^^). (2.13) 

The polynomial algebra is thus 

[A,h]=h, [A,l2] = lQ^'h, (2.14) 
[h, h] = SXiQi^iH + ^A) - X^)Q{\{H + \a))S{\{H - \a) + A,) 

-Q{\{H + \a) + 2X.^)Q{\{H + \a) + X.)S{\{H - \a)). 

The order of this polynomial algebra depends of the order of polynomials Q 
and S. Polynomial algebra given by Eq.(2.11) amd Eq.(2.14) have a similar 
structure as the cubic algebra given by Eq.(1.9). 



3 Applications 

There are two Smorodinsky-Winternitz potentials [6] that allow sepa- 
ration of variables of the Schrodinger equation in Cartesian coordinates : 
V{x, y) = ^{x^ + y^) + ^ + j, and V{x, y) = ^ {Ax^ + y^) + bx+^. They are 
well known quadratically superintegrable systems and we apply the construc- 
tion to these two systems. They have in the y axis the following creation 
operators 

-'- ^ ^ '^^ 2 ^ + ^ 1) (3 1) 

It was first obtained in Ref. 6 and reobtained in a systematic study of sys- 
tems with second order ladder operators [37]. The three other creation and 
annihilation operators have the same form. The polynomial Q{Hx) is given 

by 
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The polynomial S{Hy) is also given by the Eq.(3.2) (by replacing by Hy 
and b by c). We can form with Q{H^), S{Hy), Eq.(2.10) and (2.11) the in- 
tegrals and the polynomial algebra. In the two cases the polynomial algebra 
is a cubic algebra where the generators are second, third and fourth order 
operators. Hovewer, we can form a simpler algebraic structure for these Ha- 
miltonians, a quadratic algebra where the generators are two second order 
and one third order operator [16]. The integrals obtained from the construc- 
tion of Section 2 are not necessarily the integral of the lowest possible order. 

3.1 Systems with a third order integral 

We considered well known quadratically superintegrable systems. We will 
apply the method to Potential 1,5 6 and present their integrals and polyno- 
mial algebras. The first system that we consider is the following [27,30] 

H=\Pl + \Pi + ^'(^^^ + T^2 + T^^) ■ (3-3) 
2 2 ^ 8a^ {x - a)^ {x + a)^ 

The creation operators are given by 

, d 1 ,1 ^ ^ 2 d d 1 ,1 Inx 

Aa^ dx 2a? x-a x + a dx dx 2a^ x-a x + a 

(3.4) 

We construct the known integral B of order 3 [27,30] from Eq.(2.10) {B = 
I[ = ^^-^/i). We have A = The deformed oscillator algebras in the x 
and y axis are given by Eq.(2.7) and (2.8) with the following expressions 



S{Hy) = 2Hy 



2^' 



We get from the Eq.(2.11) of the previous section and Eq.(3.6) the following 
cubic algebra that coincide with the one found in Ref. 30. 

[AI[]=I^, [A,I^] = —I[, (3.7) 
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a? 

fe4 t6 fc6 fc8 

-S%H' + 2^A-2^H-6^. 

a2 a4 

The energy spectrum was calculated from the Fock type unitary representa- 
tions [30]. In Section 4, we will extend this algebraic method of calculating 
the energy spectrum of superintegrable systems with a polynomial algebra 
of order seven. 



3.2 Potential 6 

The next system that we consider is an Hamiltonian for which no polyno- 
mial algebra were found from the second and third order integrals of motion. 
We will show how the procedure of Section 2 will allow us to find a quintic 
algebra. The Hamiltonian 

(3.8) 

has the following second order integral A = — Hy and third order integral 

i? = 2L^-3a^({L,P,^} + {L,P^^}) + ^{(12%+3(^)(a;^V)+24y^^^ 

(3.9) 

lUyx^ , (3^-jy)ix^ + a^ jy^ - x^)iy^ + 

73 + 48?/ 7-T , iirj 



— (x^ — a2)2 (y2 _ q2^ 



- — {(124x + 3(— )(?/2 + x^) + 24x^^ ^ ^ 



4 '"^ ' ^a^ (x^ — a^) y'^ — cP- 

i^i^-^)^y^ (x^-y^)(x^ + a^ ^ 



The creation operators are given by Eq.(3.4) in the x and y axis (by replacing 
X by y). The polynomial algebras in the x and y axis are given by Eq.(2.7) 
and (2.8) with Q{Hx) and S{Uy) given by Eq.(3.6) (by replacing by Hy 
for S{Hy)). We have A = The integrals of motion are given by Eq.(2.10) 
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{I[ = ^x^/i). Thus, we obtain with the Eq.(2.11) 



[A,I[] = Ii [A,I^] = —I[ (3.10) 



16 2 

H rAH^ -A + -^AH. 

8a4 2a4 IGa^ 



The integrals I[ is related to integrals A, B and C by 

-1 3/7^ 

3.3 Potential 5 

The Hamiltonian 

2 2 ^ 8a^ (a; + a)2 (x - a)^ 

has a quadratic A = — Hy and a cubic integrals 

i.^2L3-3a^{L,Pj + .^{^-^g^ (3.13) 

, (x^-3a^) (3j/^-8a^) 2y^ , V(x^ + g^) 
12a2 x2-a2+ (x^ - a2)2 

The integral A,B and their commutator C do not generate a cubic algebra. 
We will construct other integrals of motion from the creation and annihilation 
operators. The creation operators are given by Eq.(3.1)(by replacing x by y 
and ui by a) and (3.4). We have in the x and y axis polynomial algebras given 
by Eq.(2.7) and (2.8) with A = ^ and 

g(i..)^2^3_T^^.^r^^^^l5^, (3.14) 

5™ = ^^?-^^.-^ • 
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We obtain with the Eq.(2.14) and (3.31) the following polynomial algebra 
with integrals given by Eq.(2.10)( I[ = a^h) 

[A,I[]=I'„ [A/^] = ^/(, (3.15) 
^ ^' 2^ ~ 64aio 64^8 IQa^ ^ IQa^ 4.0? ^ 4 



64 ' 64 64 ' ' 16 16 64 ' 

+Al r- \ a^H^h^^ )+A?{ tH o?H^h^ 

^Ua? 64 16 16 ' ^Ua^ 8a^ 8 2 

^64a8 16a6 16a4 2a2 4 ^ 

The integral /( is of order 7 and the integral /2 is of order 8. 



4 Realizations of polynomial algebras 

In the previous Section, we generated polynomial algebras of many sys- 
tems. These algebras were cubic, quintic and seventh order algebras. In earlier 
articles it was demonstrated that quadratic [16] and cubic [30] algebras can 
be realized as deformed oscillator algebras [38]. Deformed oscillator algebras 
allow to construct Fock type representations and obtain the energy spec- 
trum. We will show that we can construct similar realizations. We consider 
the following polynomial algebra of order seven 

[A,B] = C, [A,C] = 6B, 

[B,C]=mA'^ + nA^ + fiA^ + uA^ + aA^ + PA^ + -fA + e . (4.1) 

where A and B are integrals and thus commute with the Hamiltonian H. 
The structure constants m, n, i/, a, (3, 7 and e are polynomials of the Ha- 
miltonian. We do not impose a order to these polynomials and only make the 
hypothesis that integrals generate an algebra of the form given by Eq.(4.1). 
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The Casimir operator satisfies 

[K,A] = [K,B] = [K,C] = 0, (4.2) 

and this implies 

K = C^-6B^ + + '^nA' + (| + ^^m)^^ + (^z. + 6n)A'^ (4.3) 

+ (^ + _ )^4 (2 2 _ 1^2 

^2 6^ 12 ^ ^S'^ 3 3 ^ 

+ (^ _ 1^2 153^)^2 + -6f3 + — - —)A. 

^2 6 3^ 21 15^ 

The order of the Casimir operator depends of the order of A and B. Ultima- 
tely, the Casimir operator is written in terms of the Hamiltonian. There is a 
realization in terms of deformed oscillator algebras of the form 

A = 5{N + u), B = b^ + b, (4.4) 
where u is an arbitrary constant. Where {A^, 6, 6^} satisfy 

[N,b] = -b, [N,b^=b\ bb^ = <^{N+l), b^b = <!>{N). (4.5) 



With the third relation of the seventh order algebra given by Eq.(4.1) and 
the Casimir operator given by Eq.(4.3) we find 
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$(iV) = --S-^{N + ur+{^-'^){N + uy + {^ + -m6-'-'-^){N + u)' 

(4.6) 



,3 



.uS^ ii5'^ 1 ,.5, ,^ ,a5 5ii6'^ uSi 7 , 

+ (^ - V + 4"'^'^^'^ + ^) + + ^ - — - 48'^ "^^^'^ + ^) 

H^+'^-^-^Sln)iN+urH--^-^ + ^-^ + -S'm)iN+ur 
^ 6 6 4 12 ^ ; g 24 4 4 24 ' 

+ ( — r H - (52z/)(A^ + u) ^ -. 

^25i 12 4 84 60 ^ 4(5^ 4(5 

To obtain unitary representations we should impose three constraints on the 

structure function 

$(p + l,Mi, A;) = 0, $(0,11, A;) = 0, ^(x) > 0, V a; > . (4.7) 
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5 Potential 5 and 6 



The Eq.(4.6) gives the structure function in terms of the parafermio- 
nic number and the structure constants. The three conditions given by 
Eq.(4.7) provide a method to obtain the energy spectrum. From the results 
of Section 4 and 3, we can found the unitary representations and the corres- 
ponding energy spectrum for the Potential 5 and 6. 

5.1 Potential 6 and quintic algebras 

The algebra of Potential 6 is given by the Eq.(3.10) is a particular case 
of the one given by the Eq.(4.1). The structure constant are 

- /i=-T7r, iy = f3 = e = 0, (5.1) 



a 



4 



16 ' 



2 8a* a"^ 2a* a° loa* 

This quintic algebra is generated by integrals A, I[ and I2 respectively of 
order 2, 5 and 6. We can write the Casimir operator given by Eq.(4.3) as a 
polynomial of the Hamiltonian only 

K = -Ah'H' + - ^-\h' - + -^H' + 4rr^ - 

Aa^ a^^ 16a^^ 

(5.2) 

We can found with the Eq.(4.6) the structure function and factorize it in the 
following way 

((^ + -)-(^-^))((^ + «)-(^ + ^)) 

{{x + u)-{^ + ^))((x + u)- (^^ + ^)). 

To obtain unitary representations we should impose three constraints given 
by Eq.(4.7). There are four solutions for a = iuQ. G M. Let us present these 
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unitary representations with the corresponding constant u and energy spec- 
trum : 



Case with m = ^ + ^ 



^0 

$i(a;) = — gx(x + 2)(x + 3)(p + 4 - + 3 - + 1 - x). (5.5 
Case with U2 = ^ — | 

^six) = -r^x(x - 2)(p + 4 - + 3 - 1 - x),p = 0, 1. (5.7 

Case with U3 - ^ 



ft2 2 

2a; 



^3 = (5.8 



^0 



4ag 



$3(x) = — x(x - 3)(p+ 1 - x)(p + 3 - 2;)(p + 4 - = 0, 1,2. (5.9 

;i2(p-3) 
4 



- (5.10 



$4(x) = — gx(x — 3)(p + 1 — x)(p — 2 — x)(p — = 0, (5.11 
4aQ 

We must also exclude spurious states with other conditions. One of them 
consists in £^ > min V . We have unitary representations valid only for p=0, 
p=0,l and p=0,l,2. Such solutions were also found in the context of cubic 
algebras. This phenomenon is related to zero modes, singlet state, doublet 
states and higher order supersymmetric quantum mechanics [30,31]. A singlet 
state is annihilated by the annhilation and creation operators. The energy 
spectrum is confirmed by the results obtained from supersymmetric quantum 
mechanics. 

There is one solution for the case a G M with u = ~^2^ + §■ The unitary 
representation is 

= P>3, (6.12) 

$i(a;) = — -x{p + I - x){x + ?>){p + A - x){p + 2 - x). (5.13) 
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5.2 Potential 5 and polynomial algebras of order 7 

The algebra of the Potential 5 is given by (3.15). This case belongs to the 
one given by Eq.(4.1). The structure constants are obtained by comparing 
the Eq.(3.15) and (4.1). 

This seventh order algebra is generated by integrals A, I[ and I2 respectively 
of order 2, 7 and 8. The Casimir operator given by Eq.(4.3) can be written 
as a function of the Hamiltonian only 

K = a^H' - Aa'ffH' + m'^H^ + ^H'' - ^^h' (5.14) 
2Qin}^ . 133;ii2 275^11^^ 1425^1^6 

H 7. H — 7^n TTtH H TTT ■ 

4a6 16a8 16aio 256ai2 

The Eq.(4.6) give us the structure function and we can factorize it in the 
following way 

^, , ,4/i^2 \ CL^E.., , 1 ci^E.. , , 

$(^) = {^){x^u- (-- -^)){x^u- (-- + ^)) (5.15) 

, ,1 a^E... ,3 a^E... ,5 c?E.., ,5 a^E.. 

, ,5 a^E.., ,7 a^E.. 

Let us present the solutions for the case a = ao«, ao ^ 
There is two solutions for u = | + 

El = ^7-^, $i(x) = (f^)x(p+l-x)(2p+3-2x)(2p+5-2x)2 (5.16) 
2ao 4aQ 

(p + 2 - x)(p + 3 - x)(3 + 2x). 

(p + 1 - x)(p + 2 - x)(p + 3 - 2:)(-3 + 2x) 

The representation $2(2;) is valid for p=0. We confirmed these energy levels 
with the results obtained using the separability in Cartesian coordinates and 
the SUSYQM. 



^2 = ^^^^^X^' '^'2(2^) = (^)(3 + 2p-2x)(5 + 2p-2x)2 (5.17) 
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There is one solution for the case a G M. For u = 2 — 



we get 




x) 



(5.18) 



(5.19) 



6 Conclusion 

In this article, we have showed how we can construct integrals of mo- 
tion for two-dimensional Hamiltonians of the form given by Eq.(2.1) from 
the creation and annihilation operators of one-dimensional Hamiltonians 
and Hy. We construct from these integrals of motion higher order polynomial 
algebras. We presented the algebras for two cases ! Xx — Ay and — Xy. 
These polynomial algebra are given by the Eq.(2.11) and(2.14). These re- 
sults, we explained the particular form of the cubic algebra obtained in Ref. 
30 and given by Eq.(1.6). 

One important result of this article is the construction of higher polyno- 
mial algebra for the Potential 5 and 6. They are respectively seventh order 
and quintic algebras. The quantum case appears to be more richer than the 
classical one in terms of algebraic structures. The integrals of motion of all 
classical superintegrable systems with a second and a third order integral in 
E2 generate cubic Poisson algebra [28]. 

We have also studied the realization in terms of deformed oscillator algebras 
of a class of polynomial algebras of the seventh order. These results allowed 
us to obtain the structure function for the Potential 5 and 6 and to obtain 
unitary representations with their corresponding energy spectrum. We co- 
robored these results with those obtained by the means of supersymmetric 
quantum mechanics [30]. 

These results can be generalized in higher dimensions. The method of the 
Section 2 can also be used to generate new superintegrable systems from 
known one-dimensional systems. 

Acknowledgments The research of I.M. was supported by a postdoctoral 
research fellowship from FQRNT of Quebec. The author thanks P.Winternitz 
for very helpful comments and discussions. 



15 



7 References 



1. V.Fock, Z.Phys. 98, 145-154 (1935). 

2. V.Bargmann, Z.Phys. 99, 576-582 (1936). 

3. J.M.Jauch and E.L.Hill, Phys.Rev. 57, 641-645 (1940). 

4. M.Moshinsky and Yu.F.Smirnov, The Harmonic Oscillator In Modern Phy- 
sics, (Harwood, Amsterdam, 1966). 

5. J.Fris, V.Mandrosov, Ya.A.Smorodinsky, M.Uhlir and P.Winternitz, Phys.Lett. 
16, 354-356 (1965). 

6. P.Winternitz, Ya.A.Smorodinsky, M.Uhlir and I.Fris, Yad.Fiz. 4, 625-635 
(1966). (English translation in Sov. J.Nucl.Phys. 4, 444-450 (1967)). 7. A.Makarov, 
Kh. Valiev, Ya.A.Smorodinsky and P.Winternitz, Nuovo Cim. A52, 1061- 
1084 (1967). 

8. N.W.Evans, Phys.Rev. A41, 5666-5676 (1990), J.Math.Phys. 32, 3369-3375 
(1991). 

9. E.G.Kalnins, J.M.Kress, W.Miller Jr and P.Winternitz, J.Math.Phys. 44 
(12) 5811-5848 (2003). 

10. E.G.Kalnins, W.Miller Jr and G.S.Pogosyan, J.Math.Phys. A34, 4705- 
4720 (2001). 

11. E.G.Kalnins, J.M.Kress and W.Miller Jr, J.Math.Phys. 46, 053509 (2005),46, 
053510 (2005), 46, 103507 (2005), 47, 043514 (2006), 47, 043514 (2006). 

12. E.G.Kalnins, W.Miller Jr and G.S.Pogosyan, J.Math.Phys. 47, 033502.1- 
30 (2006), 48, 023503.1-20 (2007). 

13. J.Daboul, P.Slodowy et C.Daboul, Phys.Lett. B317, 321-328 (1993) 

14. Ya.I.Granovsky, A.S.Zhedanov et I.M.Lutzenko, Ann. Phys. (New York) 
217, 1, 1-20 (1992). 

15. P.Letourneau et L.Vinet, Ann. Phys. (New York) 243, 144-168 (1995). 

16. C.Daskaloyannis, J. Math. Phys. 42, 1100-1119 (2001). 

17. E.Witten, Nucl.Phys. B188, 513 (1981) ; E.Witten, Nucl.Phys. B202, 253- 
316 (1982). 

18. L.Gendenshtein, JETP Lett., 38, 356 (1983). 

19. B.Mielnik, J. Math. Phys. 25, 12 (1984). 

20. G. Junker, Supersymmetric Methods in Quantum and Statistical Physics, 
Springer, New York, (1995). 

21. M.S.Plyushchay, Annals Phys. (N.Y.) 245, 339 (1996) ; M.Plyushchay, 
Int. J.Mod.Phys. A15, 3679 (2000) ; F.Correa and M.Plyushchay, Annals Phys. 
322, 2493 (2007). 

22. F.Cooper, A.Khare et U.Sukhatme, Phys.Rept. 251, 267-385 (1995). 



16 



23. A.Khare and R.K.Bhaduri, Am.J.Phys 62, 1008-1014 (1994) 

24. I.F.Marquez, J.Negro and L.M.Nieto, J.Phys. A :Math. Gen. 31, 4115 
(1998). 

25. D.J.Fernandez C,V.Hussin and L.M.Nieto, J.Phys.A :Math. Gen. 27, 3547 
(1994) ; D.J. Fernandez, V.Hussin and O.Rosas-Ortiz, J.Phys.A : Math.Theor. 
40 6491-6511 (2007) 

26. S. Gravel and P.Winternitz, J.Math.Phys. 43 (12), 5902 (2002). 

27. S. Gravel, J.Math.Phys. 45 (3), 1003-1019 (2004). 

28. I.Marquette and P.Winternitz, J.Math.Phys. 48 (1) 012902 (2007). 

29. I.Marquette and P.Winternitz, J. Phys. A : Math. Theor. 41, 304031 
(2008). 

30. I.Marquette, J. Math. Phys. 50, 012101 (2009). 

31. I.Marquette, J.Math.Phys. 50 095202 (2009). 

32. J.Schwinger, AEG Report NYO-3071 (1952). 

33. V.A.Dulock and H.V.Mcintosh, Am. J.Phys. 33, 109 (1965). 

34. A.Cisneros an d H.V.Mcintosh, J.Math.P hys. 11 3 870 (1970). 

35. V.Sunilkumar, |arXiv :math-ph/ 02030"47 (2002). 

36. D.Bonatsos, C.Daskaloyannis, Prog.Part.Nucl.Phys. 43, 537 (1999). 

37. C.P.Boyer and W.Miller Jr., J.Math.Phys. 15, 9 (1974). 

38. C.Daskaloyannis, J.Phys.A : Math.Gen 24, L789-L794 (1991). 



17 



